RADC-TR-90-31 8 

Interim  Technical  Report 
December  1990 


AD-A232  039 


ADVANCED  MODELS  FOR 
CLUTTER 


Virginia  Polytechnic  Institute  and  State  University 
Dr.  Gary  Brown 


APPROVED  FOR  PUBLIC  RELEASE;  DISTRIBUTION  UNLIMITED. 


Rome  Air  Development  Center 
Air  Force  Systems  Command 
Griffiss  Air  Force  Bap~  K,v  13441-5700 


This  report  has  been  reviewed  by  the  RADC  Public  Affairs  Division  (PA) 
and  is  releasable  to  the  National  Technical  Information  Service  (NTIS) .  At 
NTIS  it  will  be  releasable  to  the  general  public,  including  foreign  nations* 

RADC-TR-90-318  has  been  reviewed  and  is  approved  for  publication. 


ROBERT  J.  PAPA 
Project  Engineer 


APPROVED : 


Director  of  Electromagnetics 


FOR  THE  COMMANDER: 


Directorate  of  Plans  &  Programs 


If  your  address  has  changed  or  if  you  wish  to  be  removed  from  the  RADC 
mailing  list,  or  if  the  addressee  is  no  longer  employed  by  your  organization, 
please  notify  RADC  ( EECE  )  Hanscom  AFB  MA  01731-5000.  This  will  assist  us  in 
maintaining  a  current  mailing  list. 

Do  not  return  copies  of  this  report  unless  contractual  obligations  or  notices 
on  a  specific  document  require  that  it  be  returned. 


REPORT  DOCUMENTATION  PAGE 


orm  Approved 
OMB  NO.  0704-0188 


Piifc  rapotrgbiKkn fa  trfccoMknt*  narration  katirad  to ami(pi  ftou  pa  raoorak  mxtiigtrk  ttm  fa  n.»«-q  rmnxarm  i— rchnrj  cogng  c Mt»  taicw 
LI"-'  jirrt— tirTTlrtittk  iMTwikit  rf  ifimaii  Sana  cajtttwii  iigrtiq  tfw  omkn  momma  a  try  arm  mtota  of  tna 

ccOcoon  a  rtarrMxTX  rxjxtiTQ  njggMtigna  fa  gductitg  ffq  Oman,  to  tti^kT^on  H— jgjaari  Taikm  OfremafarfantionOpoaoni  anoftaparts.  1215  Jtif»aor> 
OwH^nif,  s.a-»au  Aitiigon  ***»-«"»  «-rt tr tfi  "»n  n*  Mo  ajtiiat  «ndBudgfc  Pioaxnhi  HaOuaBn  Ptqaa  (C7Q«^»,w»angor\  0C20sm _ 


1 .  AGENCY  USE  ONLY  (Laava  Blank) 


2  REPORT  DATE 

December  1990 


a  REPORT  TYPE  AND  DATES  COVERED 
Interim  Jan  89  to  Feb  90 


4.  TITLE  AND  SUBTITLE 

ADVANCED  MODELS  FOR  CLUTTER 

5.  FUNDING  NUMBERS 

C  -  "19628-S9-K-C010 

PE  -  61102F 

&  AUTHOR(S) 

Dr.  Gary  Brown 

PR  -  23C5 

TA  -  J4 

WU  -  64 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS (ES) 

Virginia  Polytechnic  Institute  and  State  University 

Bradley  Department  of  Electrical  Engineering 

Blacksburg  VA  24061 

a  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

N/A 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADORESS(E» 

Rome  Air  Development  Center  (EECE) 

Hanscom  AFB  MA  01731-5000 

10.  SPONSORING/MONITORING 

AGENCY  REPORT  NUMBER 

RADC-TR-90-318 

1 1 .  SUPPLEMENTARY  NOTES 

RADC  Project  Engineer:  Robert  J.  Papa/EECE/{617)  377-3585 

1 2a.  D&TRBUTION/AVAIABUTY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

12b.  DISTRIBUTION  CODE 

ia  ABSTRACT («*o*runaoo»*»«0) 

This  report  summarized  three  Important  results  obtained  during  the  reporting  period. 
First,  it  is  shown  that  the  Luneburg-Kline  asymptotic  series  for  the  current  induced 
on  a  conducting  rough  surface  always  produces  an  error  in  the  shadow  zone.  That  is, 
the  L-K  series  yields  zero  net  current  in  the  shadow  zone  for  all  terms  in  the 
asymptotic  series.  This  shortcoming  limits  the  usefulness  of  the  L-K  series  for 
recovering  diffraction  effects.  The  method  of  smoothing  is  used  to  determine  the 
fluctuating  component  of  the  field  scattered  by  a  randomly  positioned  collection  of 
discrete  scatterers  to  the  level  of  the  distorted  wave  Born  approximation.  Of 
particular  Importance  here  is  the  role  of  the  coherent  or  average  field  in  the  medium 
as  a  source  for  the  fluctuating  field.  Finally,  it  is  shown  how  to  go  about  develop¬ 
ing  an  improved  high  frequency  approximation  for  scattering  from  rough  surfaces.  In 
particular,  it  is  shown  that  the  shadowed  Kirchhoff  approximation  should  comprise  the 
iterate  in  the  basic  integral  equation  for  the  current  because  it  contains  an  infinite 
resuomatlon  of  all  of  the  shadowing  effects. 


14.  SUBJECT  TERMS 
UHF  Clutter 
Clutter 


Terrain 

Statistics 


Scattering 


tapwcecooe 


7.  SECURITY  CLASSHCATION 
Of  REPORT 
UNCLASSIFIED 


1  a  SECURITY  CI_AS$HCATION  19  SECURTTY  CLASSHCATTON  20.  LIMITATION  OF  ABSTRACT 
OF  THM  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  UL 


TABLE  OF  CONTENTS 


Page 

1.  Introduction  1 

2.  An  Inherent  Limitation  of  the  2 

Luneburg— Kline  Representation  for  the 

Current  on  a  Conducting  Body 


2.1 

Background 

3 

2.2 

Analysis 

4 

2.3 

Discussion 

11 

2.4 

References 

14 

An  Application  of  Smoothing  to  Wave 

Propagation  Through  Discrete  Random  Media 

15 

3.1 

Background 

15 

3.2 

Analysis 

16 

3.3 

Discussion 

21 

3.4 

Summary 

24 

3.5 

References 

25 

A  Distorted  Wave  Born  Approximation  for  High 
Frequency  Scattering  from  Rough  Surfaces 

27 

4.1 

Background 

27 

4.2 

Analysis 

28 

4.2.1  The  Distorted  Wave  Born 

30 

Approximation  (DWBA) 

4.2.2  Application  to  High 

32 

Frequency  Scattering 

u\s) 

4.3 

Conclusions 

,/ 

33 

4.4 

References 

34 

Summary 

35 

i 


Accession  For 


NTIS  GRAAI  Sr 
DTIC  TAB  □ 
Unannounced  □ 
Justification 


By - 

Distribution/ 


Diet 


Availability  Codes 
i A vail  and/or 
Special 


.1 


1. 


Introduction 


The  primary  intent  of  this  research  is  to  develop  new  models  for  estimating  clutter 
levels  from  natural  terrain  with  particular  application  to  surveillance  radar.  Specific 
attention  under  this  contract  is  to  be  directed  toward  estimating  the  effects  of  subsurface 
volume  scattering.  For  radar  frequencies  above  s— band,  most  natural  surface  materials  do 
not  permit  much  penetration  into  the  subsurface  volume.  However,  for  frequencies  in  the 
UHF  range,  there  is  the  possibility  that  some  penetration  of  the  incident  electromagnetic 
energy  into  the  subsurface  volume  will  be  possible.  A  good  example  of  this  was  provided 
by  NASA's  Shuttle  Imaging  Radar  (SIR— B)  mission  which  detected  ancient  stream  beds 
beneath  the  desert  in  Africa. 

The  penetration  of  the  radar  signal  into  the  volume  beneath  the  surface  is  important 
because  of  the  potential  scattering  of  this  energy  by  discontinuities  in  the  dielectric 
constant  of  the  subsurface  material.  These  abrupt  changes  could  be  caused  by  rocks, 
moisture,  pockets,  tunnels,  and  vegetation  roots.  If  these  materials  scatter  energy  back  up 
through  the  surface,  they  comprise  yet  another  source  of  clutter  which  must  be  accounted 
for.  This  is  the  motivation  for  considering  the  possibility  of  subsurface  scattering  and 
attempting  to  model  it. 

The  second  chapter  presents  the  results  of  our  efforts  to  use  a  Lrneburg— Kline 
representation  to  improve  on  the  physical  optics  approximation  for  the  current  induced  on 
a  rough  surface.  We  find  the  L— 1\  representation  to  be  flawed  because  of  its  failure  to 
properly  represent  the  current  in  the  optical  shadow  regions  of  the  rough  surface.  In 
particular,  the  L— 1<  representation  predicts  that  the  current  will  be  exactly  zero  in  the 
shadow  regions  irrespective  of  frequency. 

The  third  chapter  discusses  a  distorted  wave  Born  approximation  for  the  fluctuating 
field  in  a  discrete  random  media.  This  is  the  approximation  that  we  will  be  using  to  model 
the  scattering  by  subsurface  objects,  so  it  is  very  important  to  this  research. 
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The  fourth  chapter  discusses  a  way  to  improve  on  the  standard  physical  optics 
approximation  for  the  field  scattered  by  a  rough  surface  in  the  high  frequency  limit.  This 
is  the  method  that  should  be  used  rather  than  the  Luneburg— Kline  representation. 
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2.  An  Inherent  Limitation  of  the  Luncburg— Kline  Representation  for  the  Current  on  a 

Conducting  Body 

2.1  Background 

In  the  field  of  scattering  from  randomly  rough  surfaces,  there  is  a  very  pressing  need 
to  develop  tractable  high  frequency  approximations  which  are  an  improvement  on  physical 
optics.  Of  course,  the  first  candidate  that  comes  to  mind  is  the  Fock  theory  of  shadow 
boundary  diffraction  [2.1].  However,  because  of  its  dependence  upon  certain  canonical 
geometries,  it  is  not  obvious  how  the  Fock  theory  could  be  beneficially  applied  to  the 
randomly  rough  surface  problem.  This  same  statement  applies  to  most  of  the  other 
theories  that  provide  some  improvement  on  physical  optics  [2.2]. 

Recently,  Lee  [2.3]  corrected  and  extended  the  earlier  work  of  Schensted  [2.4]  in 
applying  a  Luneburg— Kline  formalism  to  scattering  from  a  perfectly  conducting  body. 
Ansorge  [2.5]  used  the  same  methods  developed  by  Lee  to  obtain  results  for  scattering  from 
dielectric  bodies.  What  Lee  was  able  to  do  was  to  obtain  a  kQ  correction  to  the  physical 
optics  approximation  for  the  current  induced  on  a  peiiectlv  conducting  body  by  an  incident 
plane  wave.  Although  the  results  were  very  complicated  from  an  algebraic  point  of  view, 
they  did  not  require  any  special  canonical  scatterer  geometry.  Lee  obtained  these  results 
by  satisfying  boundary  conditions  on  the  surface  for  each  field  vector  expansion  coefficient 
in  the  Luneburg— Kline  series  of  inverse  powers  of  wavenumber  (kQ).  Lee  noted  that  the 
derivation  of  higher  order  correction  terms  was  possible  but  algebraically  quite  involved. 

Lee's  results  have  clearly  shown  that  the  Luneburg— Kline  formalism  is  capable  of 
providing  wavenumber  dependent  corrections  to  the  physical  optics  current  induced  on  a 
conducting  scatterer.  The  next  most  important  question  that  needs  to  be  addressed 
concerns  the  limitations  of  this  approach  and  its  subsequent  results.  That  is,  are  there  any 
fundamental  limitations  of  this  approach  and  what  are  the  significance  of  these  relative  to 
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the  accuracy  of  the  results  produced  by  the  technique?  The  purpose  of  this  paper  is  to 
investigate  these  questions. 

The  approach  to  be  taken  in  this  paper  will  be  somewhat  different  from  the  one 
employed  by  Lee  because  he  assumed  a  well-formed  shadow  zone  and  no  multiple 
scattering  [2.3].  These  assumptions  are  not  necessary  if  one  starts  with  an  integral 
equation  for  the  current  induced  on  the  surface  such  as  the  magnetic  field  integral  equation 
(MFIE).  Thus,  the  methodology  to  be  used  in  this  paper  is  to  expand  the  current, 
normalized  by  the  kQ-dependent  phase  factor  contained  in  the  incident  field,  in  a 
Luneburg— Kline  series  and  then  to  use  the  MFIE  to  solve  for  the  series  coefficients.  It  will 
be  shown  that  each  coefficient  satisfies  an  integral  equation  of  the  second  kind  in  which  the 
Born  term  requires  knowledge  of  all  lower  order  coefficients.  Furthermore,  these  integral 
equations  can  all  be  solved  exactly  because  of  the  particular  form  of  the  kernel.  The 
solution  shows  what  appears  to  be  a  fundamental  limitation  of  the  Luneburg-Kline 
representation  for  the  induced  surface  current,  namely,  that  all  the  expansion  coefficients 
are  identically  zero  in  the  shadow  zone.  That  is,  none  of  the  field  penetration  into  the 
shadow  zone  that  is  known  to  occur  for  kQ  finite  can  be  recovered  using  the 
Luneburg-Kline  formalism.  This  deficiency  is  traced  to  the  asymptotic  nature  of  the  series 
and  its  inability  to  represent  creeping  wave  type  currents  in  the  shadow  zone.  The 
conclusion  of  this  paper  is  that  the  Luneburg-Kline  representation  for  the  current  induced 
on  a  perfectly  conducting  body  is  incapable  of  accurately  predicting  the  current  in  the 
shadow  zone  and  is  therefore  limited  in  its  capabilities  irrespective  of  the  number  of  terms 
retained  in  the  series. 

2.2  Analysis 

The  problem  to  be  addressed  here  is  the  determination  of  the  electric  current 
density  3  induced  on  the  surface  z  =  £(r.)  by  an  incident  magnetic  field  ft1.  (The 
generalization  of  these  results  to  a  arbitrarily  shaped  closed  body  is  not  difficult.)  The 
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surface  z  =  C(rt)  separates  free  space  (z  >  ()  from  a  perfectly  conducting  medium  (z  <  (). 

The  unit  vector  L  specifies  the  direction  of  travel  of  the  incident  field.  The  electric  surface 
current  density  3g  must  satisfy  the  magnetic  field  integral  equation  (MFIE)  as  follows; 

Jg(r)  =  2n(r)xft1(r)  +  2n(r)x  /Jg(r0)  x  VQG(r  -  ro)dsQ  (2.1) 

In  (1),  n(r)  is  the  unit  normal  to  the  surface  at  the  point  r  =  r  +  ((r.  )z  and  is  given  by 

l  L 

n(r)  =  [  -  C/  -  y  +  z]/(l  +  ex  +  C“)1/2  (2.2) 

where  ^  =  d(/dx  and  C  -  dC,/ dy  are  the  x  and  y  surface  slopes  at  the  point  r  on  the 
surface.  G(r  —  rQ)  is  the  free  space  Green's  function,  i.e. 

G(r~r0)  =  exp(-jk0|r -rQ|)/4T|r  -  r0|,  (2.3) 

and  VQ  is  the  conventional  three-dimensional  gradient  evaluated  on  the  surface  (z  =  (  and 

zo  =  Notin§  t*iat  t*ie  area  integration  over  the  surface  can  be  converted  to  one  over 
the  z  =  0  plane  through 


where  drt  =  dxQdyo,  (1)  can  be  rewritten  as  follows; 

J(f)  =  +  2ff* /t(t„)»V0G(|?-f0|)df 

0 


(2.4) 
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where 


J(T)  =  Js(f)(i  +  <“  +  C°-)I/2 

(2.5) 

and 

ft  =  -  cx*  -  V  + i 

(2.6) 

If  J(r)  can  be  determined  from  (2.4),  the  scattered  field  can  be  found 

from  the  following 

integral  expression; 

fts(ft)  =  Vxp(?)G(|ft-f|)d?t 

(2.7) 

As  a  preparatory  step  to  introducing  the  Luneburg-Kline 

expansion,  the  ky 

dependence  introduced  by  the  incident  field  is  removed.  That  is,  with 

J(T)  =  t(r)  exp  (-jttj-r). 

(2.8) 

fl'(r)  =  exp  (— jlcj-r). 

(2.9) 

and  Hj  a  constant,  (2.4)  may  be  rewritten  as  follows: 

L(r)  =  2N  x  fi.  +  2N  *  f  £(r(J)  *V0G(  Jr  -  rQ| )  exp  (j^-Ar)drt 

(2.10) 

-x  ^ 

where  Ar  -  r  —  r  .  The  purpose  of  (2.8)  is  to  remove  the  known  high  frequency  behavior 
or  dependence  on  kQ  from  the  current.  That  is,  as  kQ->oo  it  is  known  that  C(r)  is 

independent  of  kQ.  The  modified  current  t(r)  is  expanded  in  a  Luneburg-Kline  series,  i.e. 
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£(?)=  S 

n=o 


(2.1!) 


where  the  vector  expansion  coefficient  jn(r),  n  =  0.  1 .  are  independent  of  the 

electromagnetic  wavenumber  kQ.  (2.11)  is  next  substituted  into  (2.10)  and  it  is  assumed 
that  term  by  term  integration  is  permissible  so  that  the  following  results; 


S  V^o”  =  *  fr  +  2$  x  S  h  n  f  j  (r  )  x  V  G  exp  (jtc.  •  Ar)dr 

n=o  n=o  o 


(2.12) 


In  order  for  (2.12)  to  be  satisfied,  the  integral  term  must  also  have  a  Luneburg-Kline 
expansion.  That  is,  if 


tn^’  V  =  /  V'o)  x  VoG  exp  $i  ‘  A^t 


(2.13) 


then  it  may  be  written  as  follows; 


where  the  t'nm  expansion  coefficients  are  independent  of  k  .  Combining  (2.13)  and  (2.14) 
yields 

f  f 

Tn0  +  ~~  +  ~rl  +  -  =  /  Jn^o)  xVoG(l  A?D  exp  (j^i*Ar)cirt  (2.15) 

ko  o  0 

so  that  the  vector  coefficients  can  be  determined  as  follows; 
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tnQ  =  limw/  Jn(rQ)  x  ?oG(|Ar|)  exp(j^-Ar)(kt 


(2.16a) 


f„l  =  1*0.  ko  (/In  “  VoG  expatj-Afjd?,  -t, 

n  0 


(2.16b) 


in— 1  f 

T  =  lim  km  fj  x  Y  G  exp(jt<.  ■  Ar)dr  —  i  — ^ 
nm  ko-»  o  J  Jn  o  i  tQ  p=0  kp 


(2.16c) 


Substituting  (2.15)  into  (2.12)  yields 


v  j  (f) k  11  =  2ft  x  \t  +  E  S  f  (r)k  1 
Jnv  '  o  l  nnr  o 

n=o  n=o  m=o 


(2.17) 


so  t hat  equating  like  powers  of  k^  yields 


J0(r)  =  2ft  x  {fi.  +  t0Q} 


(2.18a) 


"j  1  (r)  =  2ft  x{f10(f)  +  fQ1(r)} 


(2.18b) 


Jo(r)  =  2ft  *  {'t'n(r)  +  f9n(r)  +  fn9(r)} 


Substituting  from  (2.16)  into  (2.18)  yields  the  following  sequence  of  integral  equations  for 
the  vector  expansion  coefficients; 


!0(r)  =  2ft  x  ft  +  2ft  x  lim^/jjf)  «  VQG(  |  Ar  | )  exp(jft-  Ar)drt 


(2.19a) 
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JL(r)  =  x  !>"?»  (k0[ ./t0(r 0)  x  V0G(iA?|)  exp(j^  •  A?)d?t  -tQ()]} 

o  0 

+  2N  x  Um^/T^r)  x?qG(  |  Ar  |  )exp(j^  •  Ar)dr(  (2.19b) 

o  0 

•  • 

•  • 

•  • 

This  sequence  of  integral  equations  has  a  number  of  interesting  properties.  First, 
except  for  the  source  or  Born  term,  all  the  integral  equations  are  identical  in  form  in  that 
they  appear  as 

Jp(r)  =  2$  s  (r)  +  2ft  x  Mmx  J"j  (r)  x  VQG(  j  Ar  |  )exp(jtc.  •  Ar)drt 

0  0 

p  =  0,1,2,...  (2.20) 

Also,  since  the  source  term,  2,ft  x  s  (r),  in  (2.20)  depends  on  j“Q,  jp  ...,  and  jp_p  the 
integral  equations  are  recursive.  Thus,  if  *  can  be  determined  then  it  should  be  possible 
to  determine  all  higher  order  vector  coefficients.  Because  the  vector  expansion  coefficients, 
jn(r),  are  independent  of  kQ,  it  is  further  noted  that  the  source  term  s  (r)  is  determined 
almost  completely  by  a  Luneburg— Kline  series  representation  for  the  following  integral 

/  VoG(|Ar|)exp(jt  -Ar)drt  . 

For  example,  from  (2.19a)  it  follows  that 


(2.21a) 
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while  from  (2.19b) 


Sj(f)  =  lim^  kQ 


/T0(r0)  *  V0G(  I | )  exp(j!o-Ar)drt 


iir?oo/j0(^  x  VoG(lA^)  exp(j^i •  Ar)drt 

o  0 


(2.21b) 


It  should  be  noted  that  the  presence  of  jQ(ro)  under  the  integral  signs  in  (2.21b)  has  a  very 
minimal  effect  because  the  dominant  terms  are  the  ones  which  depend  upon  kQ. 

One  very  important  point  about  the  sequence  of  integral  equations  in  (2.19)  is  that 
their  solutions  are  known.  In  fact,  (2.19a)  is  essentially  the  magnetic  field  integral  equation 
in  the  high  frequency  (k^*)  limit,  so  its  solution  is  given  by  [2.6] 


2ft  x  fij  (f  not  shadowed) 

I0(r)  = 

0  ( r  in  shadow) 

There  is  also  the  possibility  of  a  multiple  scattering  contribution  to  J0(r)  from  other  points 
on  the  surface  [2.6].  Thus,  for  the  nth  vector  expansion  coefficient,  the  solution  is  (except 
for  the  contributions  of  multiple  ray  bounces  on  the  surface) 


^ 

2N  *  sfi  (r  not  shadowed) 

Jn(r)  -  (2.22) 

0  (r  in  shadow) 

The  complete  surface  current  density  is  thus  obtained  by  combining  (2.5),  (2.8),  (2.11),  and 
(2.22),  i.e. 
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!,(?)  =  (1  +  <i  +  {p‘l/5  expf-jSj  •  f )  S  ^TT  (2'23) 

b  x  y  n=o  k" 

The  reader  is  cautioned  not  to  expect  a  one-to-one  relationship  between  a  specific 
j"  coefficient  in  (2.23)  and  the  nth  order  iterate  of  the  original  integral  equation  in  (2.4). 
It  can  be  shown  that  a  direct  relationship  does  exist  for  n  =  0  and  1  but  it  breaks  down  for 
n  >  2.  This  results  is  a  consequence  of  the  fact  that  any  one  iterate  contains,  in  general,  all 
orders  of  l/kQ. 

2.3  Discussion 

One  of  the  immediate  consequences  of  (2.23)  is  that  the  current  on  any  shadowed 
portion  of  the  surface  is  identically  zero.  This  is  obviously  a  high  frequency  approximation, 
but  the  analysis  presented  above  makes  no  explicit  approximations  and,  in  fact,  appears  to 
be  exact.  Clearly,  an  exact  analysis  cannot  lead  to  an  approximate  result.  What  is 
happening  in  this  case  is  that  the  L— K  series  is  dong  the  best  job  that  an  asymptotic  series 
can  do  in  representing  the  current  in  the  shadowed  parts  of  the  surface.  The  failure  of  the 
L— K  asymptotic  series  is  linked  to  the  fact  that  the  current  in  the  shadow  zones  of  the 
surface  cannot  be  represented  by  an  asymptotic  series  of  the  L— 1\  form.  To  prove  this, 
recall  that  the  definition  of  an  asymptotic  series  such  as  (2.11)  is  that  if  represents  the 
partial  sum  of  the  first  m  +  1  terms  *  then 

lim„k™[t(ko.q-Sml  =  0  (2.24) 


m 

£  Jn(r;/k 

n=0  n 


n 

o 
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for  all  values  ofm.  In  the  limit  as  kQ-+oo,  the  current  in  the  shadow  zones  of  the  surface  is 
zero.  Thus,  from  (2.24)  with  m  =  0,  =  0  or  jQ(r)  =  0.  For  m  =  1,  (2.24)  yields 

lim  (k  L)  =  lim  (k  $,) 

k  00  v  0  '  £  ->  °°v  0  1 ' 

o  o 

but  =  3" j/k  because  =  0,  so  this  leads  to 

J!  =  li!}.  (k0£)  (2-25) 

0 

It  is  well  known  that  the  current  in  the  shadow  zone  has  the  form  of  a  creeping  wave  [2.1] 

2/3 

ta^E  C/exp{-j[k0«+exp(-jT/3)^k0/2)l^jr  k  db  (2.26) 

where  are  the  launching  amplitudes,  b  is  the  distance  measured  along  the  surface,  ^  is  a 
root  of  the  Airy  integral  [2.1],  and  k  is  the  curvature  of  the  surface  at  the  distance  b 
measured  along  the  surface.  Thus,  in  the  limit  as  kQ-*oo,  each  term  in  the  series  of  (2.26) 
exhibits  an  exponential  dependence  on  kQ.  Furthermore,  since  the  real  part  of  this 
dependence  leads  to  an  exponential  decay  with  distance,  it  is  clear  that 

Jl  =  iij3?oo  (koC)  =  0 


and,  in  fact  all  of  the  higher  j^  's  will  vanish.  Thus,  the  only  acceptable  asymptotic  series 
for  the  current  in  the  shadow  region  is  the  null  series.  The  reason  for  this  is  contained  in 
the  definition  of  the  asymptotic  series,  e.g.  (2.24).  That  is,  the  only  series  of  the  form 


m 

£  Jn(r)/k 
n=o 


n 

o 


which  can  satisfy  (2.24)  for  all  m  in  the  shadowed  regions  of  the  surface  is  the  null  series. 
Note  that  it  is  the  constraint  imposed  by  the  form  of  the  asymptotic  series  which  dictates 
the  end  result.  Thus,  other  than  the  null  series  result,  the  current  in  the  shadow  region  does 
not  have  an  asymptotic  series  representation.  This  is  an  important  result  because  it  is  the 
first  time  (to  the  author's  knowledge)  that  the  failure  of  an  L— K  series  in  the  shadow 
region  has  been  both  demonstrated  and  explained. 

While  the  L— 1<  series  does  not  lead  to  an  exact  result  for  the  surface  current  density, 
it  still  holds  the  potential  for  providing  a  tractable  improvement  to  a  pure  geometrical 
optics  solution.  This  fact  is  demonstrated  by  Ansorge's  [2.5]  calculations  for  scattering  by 
a  dielectric  sphere  using  the  ray  optics  field  approach.  The  attractiveness  of  the  current 
approach  as  developed  here  is  due  in  large  part  to  the  fact  that  the  complete  L—K 
representation  can  be  developed  entirely  from  an  integral  of  known  functions,  i.e. 


/  VoG(|Ar|)  exp(jt-Ar)drt  . 

o 


While  a  complete  L—K  series  development  for  this  integral  is  prohibitive,  it  may  be 

—2 

possible  to  obtain  the  terms  up  to  and  including  kQ  [2.7].  One  of  the  primary  advantages 

_2  0 

of  obtaining  the  kQ  and  kQ  corrections  to  the  k^  asymptotic  expansion  of  this  integral  is 
the  recovery  of  some  of  the  cross  polarizing  properties  of  a  rough  surface  in  the  high 
frequency  (but  not  optical)  limit. 

The  results  obtained  in  this  paper  raise  an  interesting  point.  Usually,  the  first 
frequency  dependent  correction  to  physical  optics  comes  from  the  non-zero  width  of  the 
transition  zone  between  the  illuminated  and  shadowed  regions  on  the  scatterer  and  the 
propagation  of  creeping  waves  into  the  shadow  zone.  If  the  L—K  series  produces  a  result 
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that  is  identically  zero  in  the  shadow  zone,  what  is  the  physical  source  of  the  k~n,  n=l,  2, 
...  terms  in  the  illuminated  zone?  It  would  appear  that  this  may  be  due  to  the  transition 
zone  encroaching  into  the  illuminated  zone  as  the  frequency  is  decreased.  This  is  the  kind 
of  question  that  needs  to  be  answered  if  there  is  to  be  a  full  understanding  of  the 
Luneburg— Kline  asymptotic  representation. 

It  should  be  noted  that  s^r)  in  (2.21b)  has  essentially  been  computed  by 
Chaloupka  and  Meckelburg  [2.8]  in  a  very  clever  application  of  integration  by  parts.  Their 
results  should  provide  the  basic  ingredients  to  study  the  question  raised  in  the  previous 
paragraph. 
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3.  An  Application  of  Smoothing  To  Wave  Propagation 

Through  Discrete  Random  Media 

3.1  Background 

The  distorted  wave  Born  approximation  (DWBA)  is  one  means  for  improving  on 
the  standard  Born  approximation  and  >t  is  particularly  useful  when  the  latter  fails  [3.8].  It 
is  based  on  splitting  the  unknown  into  the  sum  of  two  parts  in  such  a  way  that  one  part  is 
very  accurately  known  while  the  second,  and  yet  to  be  determined  part,  is  small  compared 
to  the  known  term.  This  method  had  its  origins  in  nuclear  physics  [3.8]  and  has  recently 
been  used  in  electromagnetic  wave  propagation  through  particulate  media  to  obtain  a  first 
order  estimate  of  the  incoherent  power  scattered  by  a  randomly  positioned  collection  of 
objects  [3.6,  3.9,  3.7]. 

In  the  em  propagation  problem,  the  total  field  is  split  into  the  sum  of  an  average  or 
coherent  field  and  an  incoherent  or  fluctuating  part  having  zero  mean.  For  sparsely 
populated  media,  the  average  field  is  known  to  be  the  solution  of  the  Foldy— Twersky 
integral  equation  [3.4].  Similarly,  the  fluctuating  field  in  this  limit  is  small  compared  to 
the  average  field  provided  one  does  not  progress  too  far  into  the  particulate  medium. 
Thus,  assuming  that  one  can  solve  the  Foldy-Twersky  equation  for  the  average  field  and 
that  one  restricts  interest  to  those  regions  of  parameter  space  where  the  fluctuating  field  is 
small  compared  to  the  mean  field,  it  should  be  possible  to  use  the  DWBA  to  estimate  the 
fluctuating  field  and,  hence,  the  incoherent  scattered  power. 

In  the  previous  em  propagation  studies  [3.6,  3.9,  3.7],  it  was  assumed  that  the 
DWBA  was  equivalent  to  embedding  the  particulate  scatterers  in  the  average  medium, 
illuminating  them  by  the  average  field,  and  then  letting  them  scatter  into  the  average 
medium,  (the  effective  parameters  of  the  average  medium  were  obtained  from  the 
wavenumber  of  the  average  field).  There  is  a  problem  with  this  approach;  the  average 
medium  only  has  meaning  relative  to  the  average  field.  That  is,  scattering  is  the  process 
whereby  the  average  field  is  converted  into  the  fluctuating  or  incoherent  field  and  this 
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latter  field  does  not  propagate  in  the  average  medium.  It  would  seem  then  that  the 
previous  studies  claiming  to  use  the  DWBA  have,  in  essence,  double  counted  for  the 
average  medium.  The  error  stems  from  postulating  how  the  DWBA  applies  to  the  problem 
rather  than  deriving  this  result  from  basic  scattering  relationships. 

The  purpose  of  this  paper  is  to  show  that  the  DWBA  follows  directly  from  an 
application  of  the  method  of  smoothing.  In  particular,  the  DWBA  will  be  shown  to  be 
equivalent  to  the  lowest  order  smoothing  approximation  for  the  fluctuating  field. 
Furthermore,  this  result  will  confirm  the  suspicion  that  previous  results  have  essentially 
double  counted  for  the  average  medium.  Given  the  conditions  under  which  the  DWBA  is 
applicable,  a  similarly  constrained  equation  for  the  average  field  will  also  be  obtained.  In 
the  absence  of  correlations  between  scatterers,  this  equation  reduces  to  the  Foldy-Twersky 
integral  equation.  All  of  these  results,  for  both  the  average  and  fluctuating  fields,  are 
obtained  by  a  straightforward  application  of  the  method  of  smoothing  to  the  discrete  media 
and  do  not  require  the  introduction  of  an  equivalent  continuous  random  media  [3.5]. 


3.2  Analysis 

The  approach  to  be  used  herein  is  the  application  of  the  method  of  smoothing  to  the 
integral  equation  describing  the  scattering  process  in  the  discrete  random  medium.  It 
should  be  noted  that  the  same  analytical  procedures  apply  equally  well  to  the  problems  of 
propagation  in  a  continuous  random  medium  or  scattering  by  a  random  interface.  This 
generality  of  the  procedure  follows  from  the  fact  that  smoothing  can  be  applied  to  any 
Fredholm  integral  equation  of  the  second  kind  describing  a  stochastic  process  [3.3]. 

Consider  a  volume  V  in  which  there  are  N  randomly  located  objects  having  possibly 


random  volumes  (Vn).  orientations  (Dn),  and  relative  dielectric  constants 


and  where 


n  =  1,  2,  ...,  N.  The  total  electric  field  at  any  point  in  space  (£t)  may  be  expressed  as  the 
sum  of  the  incident  field  (Ej)  which  would  exist  in  the  absence  of  the  objects  and  the 
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scattered  field  due  to  the  objects  (fig),  e.g. 


et(r)  =  fi.(r)  +  fi8(?) 


(3.1) 


The  scattered  field  may  be  related  to  the  total  field  inside  each  object  as  follows; 


fig(?)  =  iKEfiJ(rn  +  r0) 


(3-2) 


where  L  is  the  three  dimensional  integral  over  all  space,  i.e., 


«•  -  Iff  d'o 


and  the  dyadic  operator  is  given  by 


N 

E 

n=l 


5n(v  Vn’ 


G 


(3.3) 


The  i  superscript  on  denotes  the  total  field  interior  to  the  scatterer.  In  (3.3),  kQ 
is  the  free  space  wavenumber  (  —  2j/Aq)  while  Sn(?Q,  Vn,  Hn)  is  the  support  of  the  n  th 
scattering  object  whose  centroid  is  located  by  the  position  vector  r  .  Sn  is  given  by^ 


ro’  Vn’  nn 


1,  rQ  inside  Vn 
0,  rQ  outside  Vn 


(3.4) 


and  depends  only  on  the  volume  (V  )  and  orientation  (Qn)  of  the  nth  object  and  the  body 
centered  position  vector  r  .  The  dyadic  G  is  given  by 
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a  =  -  p.v.  [I  +  k~2  v0v0]  « [t  -  ?n  -  ?0]  +  I<[t  -  f „  -  ?„]  /3k2 


where  I  is  the  unit  dyad,  ^(  • )  is  the  three  dimensional  delta  distribution,  and  P.V.  denotes 
the  principle  value  which  excludes  a  small  spherical  volume  centered  at  r  —  rn  —  rQ  =  0. 
The  g  in  (3.5)  is  the  free  space  scalar  Green's  function, 


§r-rn-ro  =exP  -jkolr-rn-roM/47rlr-rn-ro 


Substituting  (3.2)  into  (3.1)  yields  the  following  integral  equation  of  the  second  kind 
for  the  total  field  in  the  presence  of  the  scatterers; 

fil(f)-fi,(h  +  Lftj;[t.tn+r0]fi|[rI1+to]  <37> 

The  first  step  in  applying  the  method  of  smoothing  to  (3.1)  is  to  decompose  into  the 
sum  of  an  average  or  mean  value  <£t>  and  a  zero  mean  fluctuating  part  .  That 
is,  with 


et  =  <et>  + 


(3.7)  becomes 


<Et>  +  a£t  =  fi.  +  LI<S  <  fij>  +  LKS 


and  averaging  this  equation  yields 
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(3.10) 


<£.>  =  fi.  +  LP  lf>  -*1 

l  1 


h  <e;> 


+  LP 


where  P  =  <  >  is  the  averaging  operator.  It  should  be  noted  that  in  (3.10), 


ks  <  e;> 


^  <  K^>  <fi|> 


because  it  is  implied  that  when  <£j>  appears  next  to  R^,  it  contains  a  dependence  on  the 
random  positions  of  the  scattering  objects,  as  (3.7)  shows  explicitly.  Subtracting  (3.10) 
from  (3.9)  gives  the  following  result; 

=  L(l-P)  Ks  <£[>  +  L(l-P)  Ks  (3.11) 

The  fluctuating  part  of  the  total  field  is  also  equal  to  the  fluctuating  part  of  the  scattered 
field  because 


et  =  +  ®s  <3-12> 

and  the  incident  field,  £.,  has  no  fluctuating  part  so  that  <*>£.  =  6$,  . 

1  L  S 

If  the  point  of  observation  r  is  successively  taken  to  be  inside  each  of  the  N 
scatterers,  (3.10)  and  (3.11)  will  result  in  2N  coupled  integral  equations  for  the  fields  <£|> 
and  inside  each  of  the  N  scatterers.  A  solution  of  these  equations  gives  the  fields  inside 
all  of  the  scatterers  and  then  the  average  and  fluctuating  fields  outside  of  the  scatterers  can 
be  computed  from  (3.10)  and  (3.11).  The  need  to  deal  with  interior  and  exterior  fields, 
even  in  a  statistical  description,  is  simply  a  consequence  of  the  boundary  value  nature  of 
this  problem  and  has  been  noted  previously  (3.2]. 
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The  goal  of  this  research  is  to  obtain  an  approximate  solution  of  (3.11)  and  then 
find  the  corresponding  equation  for  <£j>  that  is  self-consistent  with  this  level  of 
simplification.  The  most  straightforward  approximation  is  to  ignore  the  term 

L(l— in  (3.11)  because  it  appears  to  be  the  source  of  multiple  scattering  in  so  far 
as  ^  is  concerned.  We  use  the  word  appears  because,  without  this  term,  <5Et  will  indeed 
be  determined  by  <£^>  alone  but  <£’>  depends  on  (‘ifi’  through  (3.10).  To  resolve  this 
quandary,  it  is  necessary  to  substitute  for  in  (3.10)  using  (3.11),  i.e. 


<£t>  =  £j  +  LP 


+  LP 


I<eL(1-P)Ks 


<£!> 


+  LP 


L(l-P)KE«5£j 


(3.13) 


Now,  if  the  suspicious  multiple  scattering  term  L(1-P)Kj><$£’  is  ignored  in  both  (3.11)  and 
in  (3.13),  there  results 


6fit«(l-P)LRs<fij> 


(3.14) 


and 


<£t>  ss  £j  +  P 


1  +  L(1-P)KS 


<£! 


> 


(3.15) 


which  are  both  to  the  same  order  of  approximation  and  therefore  self-consistent 
relationships.  Note  that  in  (3.15)  <£t>  is  determined  independent  of  which  means 
that,  indeed,  (3.14)  is  a  single  scattering  approximation  for  insofar  a£  ^  is  concerned. 
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This  last  phrase  "in  so  far  as  is  concerned"  must  be  added  because  we  are  dealing  with 
two  field  quantities  here  and  we  must  be  precise  in  defining  exactly  what  we  mean  by  the 
term  multiple  scattering. 

It  is  clear  from  the  way  in  which  (3.11)  and  (3.13)  were  simplified  to  (3.14)  and 
(3.15)  that  the  latter  two  equations  for  and  <£^>  will  be  valid  whenever 

|(1-P)  LKE<fiJ>|  >  |(l-P)LK£afi{|  (3.16) 

This  inequality  says  that  the  fluctuating  part  of  a  one— time— scattered  (by  all  objects) 
average  interior  field  must  be  large  compared  to  a  one-time-scattered  (by  all  objects) 
interior  fluctuating  field.  It  would  seem  that  if  |<Ej>|  >  >  then  this  inequality  is 

surely  satisfied.  However,  there  may  be  instances  where  the  average  and  fluctuating  fields 
are  comparable  in  magnitude  but  (3.16)  is  still  satisfied.  The  point  of  this  discussion  is  to 
emphasize  that  it  is  (3.16)  that  must  be  satisfied  in  order  for  (3.14)  and  (3.15)  to  be  valid. 

3.3  Discussion 

Eqn.  (3.14)  represents  the  lowest  order  non— trivial  approximation  for  the 
fluctuating  field  scattered  in  and  out  of  the  random  media.  This  field  quantity  is  the 
source  of  the  incoherent  power  due  to  the  randomness  in  the  media.  According  to  (3.14) 
is  obtained  by  taking  the  fluctuating  part  of  the  field  resulting  from  N  objects 
scattering  the  interior  coherent  or  average  field  into  free  space.  This  result  reconfirms  the 
physical  intuition  that  the  fluctuating  field  is  caused  entirely  by  scattering  of  the  coherent 
field.  Contrary  to  what  is  postulated  in  some  earlier  works,  the  scattering  is  into  free  space 
and  not  into  an  average  medium.  Eqn. (3. 14s)  is  the  distorted  wave  Born  approximation 
(DVVBA)  for  the  fluctuating  field  in  a  random  media  comprising  discrete  scatterers.  If 
(3.14)  is  used  to  compute  the  incoherent  power  scattered  by  a  half-space  filled  with 
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randomly  positioned  scatterers,  it  can  be  shown  that  whenever  there  is  a  phase  retarding  or 
an  attenuation  effect  of  the  randomness  it  is  halved  relative  to  the  erroneous  double 
counting  previously  reported.  It  should  be  remembered  that  (3.14)  is  a  single  scattering 
approximation  for  the  fluctuating  field.  Hence,  the  DWBA  is  a  single  scattering 
approximation  for  the  fluctuating  field  and,  consequently,  the  incoherent  power. 
Hopefully,  this  result  answers  the  question  of  whether  the  DWBA  is  a  single  or  multiple 
scattering  approximation. 

An  interesting  sidelight  of  this  study  is  the  equation  that  was  obtained  for  the 
average  or  coherent  field,  i.e.,  eqn.  (3.15).  This  equation  was  obtained  by  retaining  the 
same  order  terms  as  were  kept  in  the  DWBA  for  This  does  not  mean  that  all  the 

terms  in  (3.15)  are  of  equal  importance.  For  example,  if  only  one  interaction  between  Kv 
and  <fi|>  is  retained,  i.e. 


<fit>  s  fij  +  PLKS  <fi|>  (3.17) 

the  resulting  equation  is  the  Foldy— Twersky  integral  equation  for  the  average  field.  To 
show  this,  it  is  first  necessary  to  split  the  averaging  operator  into  the  product  of  an  average 
over  the  random  positions  of  the  scatterers  and  the  conditional  average  over  all  other 
random  properties  of  the  objects  such  as  dielectric  constant,  volume,  and  orientation,  i.e. 

P  =  <  <  •  I  r  >  > 

" 

Thus,  (3.17)  becomes 


<£t>  =  fij  +  L  <  <  K^|rn>  <  £’> 

rn 


(3.18) 
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because  only  depends  on  the  other  random  parameters.  From  the  definition  of  K^,  i.e., 
(3.3),  the  conditional  average  becomes 


<Kklrn>  =  -ko(<Gr>_1)<S(ro’<Vn>’<fi>)>  S  G  r-r  -r 

n=l  *■ 


(3.19) 


where  <er>,  <V  >,  and  <fi>,  and  <S(-)>  represent  the  mean  relative  constant,  volume, 
orientation,  and  spatial  support,  respectively,  for  the  N  scattering  objects.  Assuming  a 
uniform  probability  density  function  for  the  r  ,  n  =  1,  2,  ...N,  i.e., 


1/V  r  inside  V 

P(rJ  = 


0  r  outside  V 
n 


(3.20) 


and  substituting  (3.19)  and  (3.20)  in  (3.18)  yields 


<Bt>  -  fii  ~  ko  (<€n>_1]  L<S>  I  ['~'n~7o]  <fit  [?n  +  ^o]>dVV 

w  n  —  I  J  ^  J 


:fit>  =  ei  -  ko  [<  V  -  ■]  '■<«> [v]  /«  [•  ~ ^ ?„] <  K  (f,  +  ?„] >  (3.2i : 


because  all  N  integrations  are  identical  and  the  series  can  be  summed.  Since 


L<S>  = 


<V  > 

n 


(3.21)  becomes 


<et> 


=&i-< 


<er>-i 


I  /<* 


n  ro 


<e: 


rn+  ro 


>  dr„dro 


(3.22) 


<Vn>  V 


where  p  =  N/V  is  the  density  of  scatterers  in  the  volume  V.  This  is  identical  to  the 
Foldy— Twersky  equation;  see  eqn.  (3.4)  of  [3.2], 


The  next  term  in  (3.15),  i.e. 


PLKSL(1-P)KS  <fij> 

clearly  involves  correlations  between  scatterers  so  it  must  be  akin  to  the  Lax-Twersky 
approximation  for  <fit>  [3.1,  3.10],  This  term  and  its  meaning  is  presently  under  study. 

3.4  Summary 

The  original  intent  of  this  work  was  to  rigorously  derive  a  low  order  approximation 
for  the  incoherent  or  fluctuating  field  in  a  discrete  random  medium  using  techniques  that 
did  not  require  the  use  of  an  equivalent  continuous  random  medium  or  were  not  largely 
numerically  oriented.  The  purpose  of  the  work  was  twofold  with  the  primary  reason  being 
to  check  some  earlier  calculations  supposedly  based  on  the  distorted  wave  Born 
approximation.  A  secondary  purpose  was  to  attempt  to  develop  a  relatively  simple,  self 
consistent,  low  order  approximation  for  both  the  fluctuating  and  average  fields.  The 
method  of  smoothing  was  chosen  to  accomplish  this  because  of  its  ability  to 
straightforwardly  provide  accurate  low  order  approximate  solutions. 

The  results  obtained  for  the  fluctuating  field,  in  the  lowest  order  of  approximation, 
were  found  to  be  equivalent  to  the  DWBA  but  without  the  double  counting  for  the  average 
medium  erroneously  postulated  in  earlier  efforts.  A  consequence  of  this  result  is  that  the 
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incoherent  power  scattered  from  a  volume  of  randomly  located  discrete  scatterers  is  larger 
than  predicted  by  these  earlier  results.  This  is  because  the  double  counting  for  the 
attenuation  effects  of  the  average  medium  effectively  masks  the  scattering  from  the  objects 
deep  in  the  medium.  A  secondary  result  of  this  study  of  the  fluctuating  or  incoherent  field 
was  that  it  could  be  clearly  shown  that  the  DWBA  does  not  contain  any  multiple 
scattering  of  the  fluctuating  field.  This  result  answers  the  old  question  of  does  the  DWBA 
include  multiple  scattering  by  first  forcing  one  to  carefully  define  multiple  scattering. 

Finally,  if  the  equation  for  the  average  or  coherent  field  is  developed  to  exactly  the 
same  order  of  approximation  as  used  in  the  fluctuating  field,  an  apparently  new  integral 
equation  is  obtained.  If  scatterer-to— scatterer  correlations  are  small  or  can  be  ignored, 
this  equation  is  shown  to  reduce  to  the  Foldy-Twersky  integral  equation.  If  correlations 
cannot  be  ignored,  the  result  appears  somewhat  similar  to  the  Lax-Twersky  equation  but 
this  situation  is  till  under  investigation. 

In  regard  to  a  secondary  purpose  of  this  paper,  it  was  found  that  the  method  of 
smoothing  is  a  very  simple  means  for  analyzing  the  propagation  characteristics  of  discrete 
random.  This  statement  applies  to  both  the  average  and  fluctuating  fields.  Furthermore, 
it  is  not  necessary  to  use  an  equivalent  continuous  media  representation  for  the  discrete 
scattering  process  because  the  discrete  character  of  the  problem  is  easily  dealt  with  by  the 
method  of  smoothing. 
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4.  A  Distorted  Wave  Bom  Approximation  for  High 

Frequency  Scattering  From  Rough  Surfaces 

4.1.  Background 

Recently,  it  has  been  shown  [4.1]  that,  for  random  surface  scattering,  first  order 
smoothing  is  capable  of  providing  an  improvement  to  the  classic  Rayleigh— Rice  first  order 
boundary  perturbation  result.  In  particular,  first  order  smoothing  yields  a  result  which 
depends  only  on  the  dominance  of  the  average  or  coherent  scattered  field  over  the  zero 
mean  fluctuating  field.  Consequently,  it  is  valid  for  small  surface  height  (relative  to  the 
electromagnetic  wavelength)  but  arbitrary  slopes,  curvatures,  etc.  The  absence  of  any 
restrictions  on  the  surface  height  derivatives  is  why  the  first  order  smoothing  result  is 
superior  to  the  boundary  perturbation  approximation;  in  the  case  of  small  slopes,  the 
former  reduces  to  the  latter  [4.1].  It  was  also  noted  in  [4.2]  that  whenever  it  is  the 
fluctuating 

scattered  field  to  be  determined,  first  order  smoothing  and  the  distorted  wave  Born 
approximation  (DWBA)  produce  equivalent  results.  This  is  true  for  scattering  by  random 
surfaces  or  by  volumes  comprising  randomly  varying  discrete  or  continuous  constitutive 
properties. 

For  random  media  problems,  first  order  smoothing  is  essentially  a  low  frequency 
approximation  because  it  is  dependent  upon  the  dominance  of  the  coherent  or  average 
propagating  or  scattered  field.  At  high  frequencies,  the  coherent  field  becomes  vanishingly 
small  so  first  order  smoothing  totally  breaks  down  in  its  accuracy.  However,  there  is  the 
possibility  of  still  using  the  distorted  wave  Born  approximation  (DWBA)  to  obtain  an 
accurate  description  of  the  scattering  process  in  this  limit.  The  key  to  success  in  using  this 
approximation  is  to  start  with  an  appropriate  Born  approximation;  it  need  not  be  the  Born 
term  appearing  in  the  integral  equation  for  the  unknown  scattered  field  or  the  current. 

The  purpose  of  this  note  is  to  apply  the  above  methodology  to  scattering  by  an 
arbitrarily  roughened  planar  conducting  surface  in  the  high  frequency  limit.  By  selecting 
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the  appropriate  Born  term  to  be  the  exact  high  frequency  scattering  limit,  which  is  known, 
it  will  be  possible  to  develop  an  approximation  which  resurrects  frequency  dependent 
diffraction  effects.  This  approximation  comes  from  the  Born  term  in  the  integral  equation 
for  the  difference  between  the  unknown  quantity  and  its  (known)  high  frequency  limiting 
behavior.  The  key  element  in  this  DWBA  methodology  is  accurate  knowledge  of  the 
behavior  of  the  scattered  field  in  the  ray  optic  limit  (kQ  =  2-k/X^  oo).  In  the  case  of 
random  surface  scattering,  the  high  frequency  problem  is  less  complicated  than  the  low 
frequency  limit  because  there  is  no  need  to  worry  about  the  average  or  coherent  scattered 
field  since  it  is  usually  vanishingly  small  as  k  ->  oo. 

The  net  result  of  this  analysis  is  a  rigorous  methodology  for  taking  known  high 
frequency  scattering  results  and  extending  them  down  to  lower  frequencies. 

4.2.  Analysis 

To  illustrate  the  method,  the  problem  of  scattering  by  an  arbitrarily  roughened 
planar  interface  will  be  considered.  Above  the  interface  defined  by  z  =  £(x,y)  is  free  space 
while  beneath  it  (z  <  0  is  a  perfectly  conducting  medium.  The  current  induced  on  the 
surface  satisfies  the  magnetic  field  integral  equation  given  by 

Js  (r)  ^(r)  ~  2n(r)  *  /VQg(r  -  rQ)  *  ?s(ro)dso  (4.1) 

where 

=  2n(?)  *  fij(f),  (4.2) 

g(r  -r0)  =  exp(— jkQ I  r  —  rQ |  )/4tt|  r  —  rQ | , 
kQ  is  the  wavenumber,  and  n(r  )  is  the  unit  normal  to  the  surface  at  the  point  r 
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(4.3) 


n(r  )  =  (-Cxx-Cyy+z)//i  +  (VC?. 


and  ftj  is  the  incident  magnetic  field.  It  is  convenient  to  convert  the  integral  in  (4.1)  to 

one  over  the  z=0  plane  and  to  multiply  both  sides  of  (4.1)  by  the  factor  exp(jk  ();  the 

resulting  current  can  then  be  Fourier  transformed  with  respect  to  x  and  y  to  yield  the  far 

zone  scattered  field.  The  factor  k  is  the  z— component  of  the  vector  wavenumber  pointing 

sz 

in  the  scattering  direction.  With  these  manipulations,  (4.1)  can  be  written  in  the  following 
operator  form; 


3(r  )  =  3’(r  )  +  LG(r,  rQ)  •  3  (rQ) 


(4.4) 


where 


J  (f  )  =  Jg  (r )  exp  (jkSz0  7l+(VC)2 

(4.5a) 

T(r)  =  (r  )exp(jkszC)v/l+(V02 

(4.5b) 

00 

L  =  ff(->dxody o 

(4.5c) 

—00 


C>  (r,  r0)  =  -2ft(r)  *  ?Qg(?  -  rQ)  exp  [jksz  (C~C)] x  (4.5d) 


and 


N(r  ) 


-v-y+z 


where  (  and  (  are  the  x  and  y-com ponents  of  surface  slope.  Eqn.  (4.4)  is  the  integral 
x  y 


29 


equation  that  will  be  considered  since  its  solution  essentially  determines  the  scattered  field. 


4.2.1  The  Distorted  Wave  Born  Approximation  (DWBA) 

The  key  step  in  developing  the  distorted  wave  Born  approximation  (DWBA)  is  to 
write  the  unknown  current  3  a 5  the  sum  of  a  known  term  (Jj,  whose  range  of  validity  is 

a. 

also  well  known,  and  an  unknown  term  (gJ)  which  is  small  in  a  region  where  ?  is  very 

a 

dominant  [4.3],  e.g. 


}  =  1  +  €? 
A 


This  sum  is  substituted  in  (4.4)  to  produce  the  following  integral  equation  for  eJ; 


E?  —  I (3*—?^)  -f  LG*^a  -f  LG*e3 


The  Born  or  source  term  in  this  equation  is 


(3  1  -  Ja)  +  kCl’^a 


substituting  this  approximation  K  (1.6)  for  g3  yields  what  is  usually  called  the  DWBA, 


J  «  J'  +  LG? 


It  is  interesting  to  note  that  (4.9)  is  potentially  very  different  from  a  first  order 
interative  solution  of  (4.4)  such  as 
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(4.10) 


}'  +  LG-}' 


This  is  because  }  and  }'  may  be  different  oi  they  may  have  differing  domains  of  validity. 

CL 

The  DVVBA  hinges  on  the  accuracy  of  }  within  a  given  region  of  parameter  space  and  it 

a 

attempts  to  extend  this  region  by  iterating  }  .  The  first  order  iteration  in  (4.10)  is 
generally  less  accurate  or  even  well  understood  because  one  seldom  knows  when  }'  by  itself 
is  accurate.  Of  course,  when  }'  xs  }  then  both  methods  produce  the  same  approximation. 

Cl 

Another  point  of  note  about  (4.9)  occurs  when  }  and  }'  differ  significantly.  In  this  case,  it 

a 

will  probably  take  many  interations  of  (4.4)  to  produce  the  same  effect  as  the  single  term 

LG-}„  in  (4.9).  Thus,  the  DWBA  is  effectively  a  resummation  of  the  standard  iterative 
solution  of  (4.4),  i.e. 


N 

E  [LG-]0}1 
n=0 


(4.11) 


where  N  may  be  infinite  or  finite. 

The  DWBA  is  not  without  its  own  set  of  drawbacks.  For  example,  having  both  an 
accurate  solution  and  knowing  the  region  of  parameter  space  over  which  the  solution  is 
valid  are  very  difficult  conditions  to  satisfy.  Further  compounding  the  problem  is  the  need 
to  know  just  how  much  (4.9)  will  extend  the  region  of  validity  of  }  .  Nevertheless,  it 
appears  to  be  a  method  for  which  a  bit  of  knowledge  can  go  a  long  way.  As  proof  of  this, 
one  need  only  turn  to  low  frequency  scattering  by  a  rough  surface  where  the  DWBA  is 
much  more  robust  than  the  boundary  perturbation  method  [4.1],  as  previously  noted. 
Finally,  although  (4.9)  was  derived  for  a  roughened  planar  surface,  it  is  clear  that  it  is 
much  more  general  than  this.  That  is,  (4.9)  is  valid  for  scattering  by  any  arbitrarily 
shaped  body  if  the  operator  L  in  (4.5c)  is  replaced  by  an  integration  over  the  x  and  y 
extent  of  the  scatterer. 
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4.2.2  Application  to  High  Frequency  Scattering 

In  the  limit  of  k  ->  oo,  it  is  well  known  [4.4,  4.5]  that  the  exact  solution  of  (4),  and 
also  (4.9),  comprises  the  sum  of  two  currents.  The  first  of  these  is  the  shadowed  Kirchhoff 
result,  e.g. 


Jshft  =  si^  (4-12) 

where  s-  is  unity  if  the  point  r  on  the  surface  is  illuminated  by  the  incident  field  and  zero 
if  it  is  shadowed  by  another  part  of  the  surface.  The  second  contribution  is  due  to  multiple 
ray  bounces  on  the  surface  and  it  will  be  denoted  as  simply  3m  because  it  has  a  rather 
complicated  form  in  general  [4.5].  It  can  be  obtained  by  solving  (4.4)  via  iteration  in  the 
limit  as  k  -<  oo  and  considering  only  the  stationary  phase  points  in  the  evaluation  of  the 
integral  term  which  give  rise  to  multiple  scattering  [4.5].  Thus,  the  appropriate  high 
frequency  asymptotic  solution  for  the  current  is  given  by 


Ja=  lim  3  =  s-f  +  3 

V00 


(4.13) 


This  result  is  the  exact  solution  of  both  (4.4)  and  (4.9)  for  the  limiting  case  of  kQ->  oo. 

Substituting  (4.13)  in  (4.9)  yields  the  following  high  frequency  DWBA  for  the 
current  induced  on  the  rough  surface 


3  —  3'  +  LG-Sj3'  +  LG  •  3m  (4.14) 

in  terms  of  the  known  limiting  forms  for  the  current.  In  this  case,  it  seems  clear  that  the 
DWBA  is  taking  the  exact  limiting  behavior  for  the  current  and  extending  it  to  lower 
frequencies.  Thus,  the  integrals  in  (4.14)  should  be  evaluated  accurately  as  possible  and 
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certainly  not  by  asymptotic  means  (as  kQ-*  oo)  for  this  would  just  lead  to  the  obvious  result 
in  (4.13).  There  is  another  interpretation  of  (4.14)  which  is  probably  more  meaningful  for 
the  rough  surface  scattering  case.  Taking  the  limit  as  k  -*  oo  can  also  be  viewed  as  being 
equivalent  to  dealing  with  rough  surfaces  which  have  no  frequency  components  that  are 
larger  than  the  electromagnetic  wavenumber  kQ.  Solving  (4.14)  by  non— asymptotic  means 
has  the  effect  of  allowing  some  frequency  components  on  the  surface  which  are  larger  than 
(the  now  finite)  kQ.  Unfortunately,  it  is  not  possible  to  predict  just  how  much  high 
frequency  (>k  )  structure  is  accurately  accounted  for  by  (4.14).  This  is  a  limitation 
common  to  most  finite  order  iteration  schemes. 

The  DWBA  is  clearly  superior  to  the  first  order  interative  Born  result  in  (4.10) 
when  dealing  with  near  grazing  incidence.  In  addition,  (4.14)  may  be  one  of  the  few 
approximate  results  capable  of  estimating  high  frequency  scattering  from  surfaces  having 
correlation  lengths  which  are  the  order  of  kQ  [4.6].  (For  these  surfaces,  the  term  high 
frequency  denotes  a  large  height  to  wavelength  ratio).  This  is  because  the  DWBA  is  based 
on  an  exact  resummation  of  multiple  scattering  effects  in  the  limit  as  kQ-t  oo.  In  contrast, 
the  integral  term  in  the  first  order  iterative  Born  approximation  in  (4.10)  may  not  even 
converge  to  a  meaningful  result  due  to  all  of  the  multiple  scatterings  taking  place  on  the 
surface. 

4.3.  Conclusions 

The  distorted  wave  Born  approximation  is  an  attractive  alternative  to  the  first 
order  iterative  Born  approximation  under  certain  conditions.  The  purpose  of  this  paper  is 
to  point  out  these  conditions  when  one  is  interested  in  high  frequency  scattering  from 
roughened  planar  conducting  surfaces.  To  accomplish  this  purpose,  a  DWBA  is  developed 
for  scattering  from  arbitrary  surfaces  and  contrasted  to  the  first  order  iterative  Born 
approximation.  In  the  case  of  high  frequency  scattering,  it  is  shown  that  the  DWBA  is 
capable  of  extending  a  ray  optic  solution  into  the  finite  frequency  domain.  The  potential 
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for  this  approximation  to  predict  scattering  from  very  rough  surfaces  is  also  of  note. 
Consequently,  it  may  well  be  capable  of  dealing  with  the  backscattering  enhancement 
problem  [4.6]  associated  with  such  surfaces. 
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5.  Summary 

The  goal  of  this  research  is  to  advance  the  understanding  of  clutter  and  to  thereby 
lead  to  better  models  for  predicting  clutter.  The  first  chapter  points  out  and  explains  a 
basic  limitation  of  the  Luneburg-Kline  expansion  for  the  current  on  any  perfectly- 
conducting  body  or  surface.  Although  not  specifically  shown,  this  result  also  applies  to 
imperfectly  conducting  surfaces.  It  is  found  that  the  Luneburg-Kline  representation  in 
inverse  powers  of  the  electromagnetic  wavelength  always  leads  to  zero  current  in  the 
optical  shadow  zone  of  the  surface.  Hence,  the  Luneburg-Kline  representation  is  limited  in 
its  ability  to  produce  the  correct  diffraction  effects  associated  with  finite  frequencies.  This 
means  that  there  is  an  uncertainty  associated  with  the  use  of  the  Luneburg-Kline 
expansion  in  that  it  is  not  clear  when  it  starts  to  produce  false  results.  Given  this 
limitation,  it  is  not  clear  that  the  L-K  representation  is  worth  further  development  or 
study. 

The  second  chapter  investigates  ways  to  most  easily  and  straightforwardly  represent 
the  scattering  by  a  volume  distribution  of  randomly  located,  shaped,  and  electrically 
constituted  objects.  The  primary  finding  of  this  chapter  is  that  the  method  of  smoothing  is 
a  very  robust  approach  for  providing  both  low  and  high  order  approximations  for  the 
average  or  mean  field  and  the  fluctuating  field.  The  average  field  is  important  because  it 
acts  as  the  source  of  the  fluctuating  field  which,  in  turn,  is  the  source  of  the  incoherent 
scattered  power.  The  smoothing  approach  shows  that,  to  lowest  order,  the  fluctuating  field 
is  formed  by  all  the  objects  in  the  medium  scatteiing  the  mean  incident  field  into  free 
space.  This  result  is  contrary  to  previous  postulated  results  wherein  the  fluctuating  field 
scatters  into  the  average  medium.  The  importance  of  this  latter  result  is  its  prediction 
that  objects  deep  in  the  scattering  volume  will  be  more  significant  contributors  to  the 
backward  scattering  process  titan  previously  thought. 

Finally,  the  last  chapter  investigates  how  one  goes  about  obtaining  an  improvement 
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to  the  Kirchhoff  approximation  for  the  surface  current.  First,  it  is  pointed  out  that  the 
pure  Kirchhoff  approximation  is  neither  a  low  frequency  or  a  high  frequency  asymptote. 
That  is,  it  lacks  the  effects  of  shadowing  to  make  it  a  true  high  frequency  result,  and  it 
does  not  contain  the  proper  polarization  dependence  for  a  low  frequency  limit. 
Consequently,  if  one  starts  with  the  Kirchhoff  approximation,  it  is  difficult,,  at  best,  to 
estimate  just  how  accurate  such  a  starting  point  is.  The  material  in  the  last  chapter 
recommends  starting  with  the  shadowed  Kirchhoff  approximation  because  it  is  an  exact 
result  in  the  optical  limit.  In  fact,  iteration  of  this  approximation  appears  to  be  one  way 
for  recovering  some  lower  frequency  diffraction  effects.  Of  course,  it  is  also  recognized  that 
the  inclusion  of  shadowing  at  the  level  of  the  current  is  easy  to  write  down,  but  it  is 
difficult  to  actually  implement.  Nevertheless,  this  chapter  does  provide  some  rationale  for 
improving  on  the  standard  Kirchhoff  approximation. 
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RADC  plans  and  executes  research,  development,  test  and 
selected  acquisition  programs  in  support  of  Command,  Control, 
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engineering  support  within  areas  of  competence  is  provided  to 
ESD  Program  Offices  (POs)  and  other  ESD  elements  to 
perform  effective  acquisition  of  CPI  systems.  The  areas  of 
technical  competence  include  communications,  command  and 
control,  battle  management  information  processing,  surveillance 
sensors,  intelligence  data  collection  and  handling,  solid  state 
sciences,  electromagnetics,  and  propagation,  and  electronic 
reliability /maintainability  and  compatibility. 


